Entanglement entropy of fermions in any dimension and the Widom conjecture 
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We show that entanglement entropy of free fermions scales faster then area law, as opposed to 
the scaling L d ~ 1 for the harmonic lattice, for example. We also suggest and provide evidence in 
support of an explicit formula for the entanglement entropy of free fermions in any dimension d, 
S ~ c(dF,dn) ■ L d ~ x logL as the size of a subsystem L — > oo, where dT is the Fermi surface and 
80. is the boundary of the region in real space. The expression for the constant c(8T, 8Q) is based 
on a conjecture due to H. Widom. We prove that a similar expression holds for the particle number 
fluctuations and use it to prove a two sided estimates on the entropy S. 



In recent years a number of parallel findings have em- 
phasized the importance of entanglement entropy [1-10]. 
Although originally studied in the context of black hole 
physics y|, this quantity also plays an important role 
in quantum information as a measure of the number of 
maximally entangled pairs that can be extracted from a 
given quantum state |2j. 

The behavior of entanglement entropy is closely re- 
lated to the criticality behavior of quantum systems: for 
a gapped system one expects an area scaling law due to a 
finite correlation length £. In the ID case |3(, this behav- 
ior changes drastically near criticality where the absence 
of gap leads to long range correlations, and so the en- 
tanglement does not saturate. Many interesting results 
have been obtained for ID models. For classes of critical 
models, where conformal field theory (CFT) methods are 
applicable, the entropy was found to exhibit a logarith- 
mic behavior, with a coefficient depending on the central 
charge of the CFT models 0, 0] , recently a modification 
of these results in case of strong disorder was found . 
The bi-partite structure of the ground state of Fermion 
models was studied in several works 0, S A con- 
nection between the entropy of spin chains and Random 
Matrix Theory was established in 0] . 

Fewer results were obtained in dimensions d > 1, al- 
though from the field theoretic point of view these are 
very interesting. Indeed, initial investigation of the en- 
tropy as a correction to the Bekenstein-Hawking entropy, 
suggested that the entropy of a scalar field is proportional 
to the boundary area for spherical or a half-space regions 
P, ^3 . Recently, it has been rigorously proved ^3] f° r a 
harmonic lattice model, that the entropy of a cube with 
side L behaves as the boundary area, i.e. as L d ~ l . 

In this Letter we examine the dependence of entan- 
glement entropy on dimension and geometry in a simple 
case of a gapless system consisting of free fermions. 

Let us summarize our main results: First, we prove 
that S ~ L d-1 logI H| for cube-like domains (© be- 
low) . We then present a heuristic argument for the more 
explicit formula S ~ |L d_1 logL. We note that results 
which are derived for cubes do not necessarily describe 
the scaling for general boundaries: indeed we find that 



S > L^ 13 for fractal-like boundaries, where (3 <E (0, 1) 
(described below) characterizes regularity of the bound- 
ary. However, the results for cube- like domains should re- 
flect the correct scaling for regions with sufficiently regu- 
lar boundaries. For general piecewise smooth boundaries 
we prove 0(L d - x logL) < S < 0{L d ~ 1 (log L) 2 ), see ®, 
@ (this estimate was independently derived in [9j for d 
dimensional cubes in the lattice case). 

Finally, making a connection with a conjecture of 
Widom we suggest an explicit geometric formula 
for the entropy as L — * oo: 
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where dT, <9f2 are the boundaries of the Fermi sea and the 
region considered, Tip j Tlx £lX(5 the unit normals to these 
boundaries. We present evidence supporting this con- 
jecture and prove a similar formula for the fluctuations 
in particle number in the subsystem, which also gives 
bounds on S. Recently the formula Q was checked nu- 
merically for 2D and 3D and a perfect agreement 
concerning both the order and the coefficient was found. 
Widom's conjecture is closely related to the problem of 
recovering data from a measurement during a finite time 
interval and in a finite frequency set. This problem, 
known as time-frequency limiting, is of basic importance 
in signal theory, and it was studied extensively 0|. It 
turns out that operators appearing in calculations of en- 
tanglement entropy for free fermions are exactly the same 
as the ones studied in which is natural since one 
studies the properties of a field in restricted sets of real 
space and momentum space. 

The ground state of a translation invariant Hamilto- 
nian describing a non-interacting fermion field (on a lat- 
tice or in the continuum), with dispersion relation e(fc), 
i s rL(fc)<c F a !|0 >• Here ep is the Fermi energy. This 
defines the Fermi sea region T = {k\e(k) < ep} in mo- 
mentum space. We also assume that the system is gap- 
less |13 • The bi-partite structure of the ground state can 
be studied by fixing a region Q in real space and com- 
puting the reduced density matrix pQ = Tr^jjd^)) p 
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FIG. 1: The Fermi sea T in momentum space, and a region 
fi in real space 



where J 7 (pi) is the fermion Fock space associated with 
the region SI (see Fig. QJ. The entanglement entropy 
S = - Tr(jr(o)) po log/on is given in this case @ by 
S(X) = Tr h(PQP) where ft = hi + ft 2 with fti(f) = 
— tlogt and /ia(t) = — (1 — t)log(l —t). Here P is a pro- 
jection operator on the modes inside the Fermi sea F, 
and Q is a projection on the region scaled by a factor 
L. The operator PQP is related to the fermion correla- 
tion function g(x — x') =< a\a x i >—< x\P\x' > so that 
<x\PQP\x' >= J Q g(x-x")g(x" -x')dx". The density 
of particles is n — vol(T), and one may rescale L appro- 
priately, as to set n = 1, which we will assume from now 
on. 

Results for cubic domains. - Consider the case of a rect- 
angular box with sides Lj i.e. SI = [0, L\] x • • • x [0, LJ\ 
and T = [0, l] d . Let Si(L) be the entropy in the ID case 
and let S be the entropy corresponding to SI, T as above. 
Then we have the following 

Theorem. Under the above assumptions 

d d 

Y d e s ^ n n ( l *) < $ < z * n ( 2 ) 

where N(Lj) is the average number of particles. Note, 
in particular, that for SI = [0,L] d 

1 / L \ d - 1 ( L \ d- 1 

-(-) Sl (L)<S<d(-) Sl (L). (3) 

Proof: Note that we can make separation of variables 
Q = ®iQj where Qj is a projection on coordinate j, and 
P factors in a similar way. Hence PQP — ®fTj where 
Tj = PjQjPj- Note the following 

Lemma. For et^ 6 [0, 1] one has 

±G(a u ■ ■ • , a d ) < h 2 {Uti a i) < G(a u ■ ■ ■ , a d ) (4) 

where G(ai,--- ,a d ) = E/=i hvfaj) Ylijtj o>i- To prove 
Q, one has to check it for two variables ai, a 2 and then 
proceed by induction. 

We observe that the eigenvalues of PQP are of the 
form a^i • ■ - ai d with a,,j being some eigenvalue of Tj. 
Writing the entropy S as M a i.i ' ' ' a i,d)i using and 



/l i(]l l =i a i) = Z)j=i ^i( a j) IIm^ a ii an d recalling that 
the average number of particles [8j AT = Tr Tj, J5J follows. 

Result shows that the entropy may be evaluated 
using the ID expressions. This reflects the compatibil- 
ity of r, SI with factorizing the fermionic modes into the 
different coordinates. It is now a matter of substituting 
the numerous results obtained in the ID case. For the 
lattice case, it follows from the many works on the sub- 
ject that for fermions on a ID lattice, or equivalently for 
an XX spin chain (via the Jordan- Wigner transforma- 
tion), Si(L) = jlogi + o(logi), see in particular [ijjj]. 
In the continuous case the same expression is obtained 
by formally substituting h(t) in the ID result of [J3 @- 

For any body composed of a union of cubes Ci 
of side Li we have, using the subadditivity of en- 
tropy [3, S({JCi) < Y^ S ( C i)> thus we have an 
upper bound that depends on the number of cubes 
needed to describe the body: using J3J we find 
SQJCi) < K^^E^log^) + o(L d ^logL M ), 
Lm = maxjLi. A lower bound proportional to 
L d ~ 1 logL m , L m = minimi follows from ©, © below. 

Scaling coefficient. - Here we derive heuristically 

S =t s (yY' 1o S L + o(L d - 1 log L) (5) 

for SI — r — [0, L] d . Since the eigenvalues of Tj 
are strictly less than one, the series Tr h%(PQP) = 
- En=2 7^T) Tr«T;) converges. By [H 
(see below for d = 1 and f(t) = t n ), as L — > oo, 
TrT " = h + ^ El I + °(log£). Hence T^sfrj*) = 
+ i) + ol^ 1 logL). Substitut- 

ing the latter in the series for ft 2 and calculating the 
sums involved we find: Tr h 2 (PQP) = f (^) d_1 logi + 
o(L d ~ 1 logi). Adding this to Tr hi(PQP) which is com- 
puted directly and gives the same value, JSJ follows [2^ |. 
Further control of the remainder term in Tr T™ as L — > 00 
is required to make this calculation rigorous [scj . 

Results for general boundaries.- We now turn to the 
case of general bounded Fermi sea V and region f2. It 
is known that the variance in particle number, given 
by (AA0 2 '= Tr PQP(1 - PQP) can be used to obtain a 
lower bound on S. 

Theorem. For general sets fl,T one has 

4(AA0 2 < S < O(logL)(AA0 2 . (6) 

We derive also an explicit formula (jSJ) which implies in 
particular that (AN) 2 = 0(L d ^ x logL). 

The proof of © in the lattice case is immediate using 
the inequalities of the form 4£(1 — t)< h(t) < e — Ct(l — 
t) loge, valid for e > 0, with C being a constant 18] . One 
substitutes the operators PQP instead of t and calculates 
the trace. Note that tr(e) ~ eL d for a finite lattice of size 
L, thus taking e < ^J^, © follows. The proof of © in 
the continuous case is new: note first that the kernel of 
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the operator PQP is given by 

< P \PQP\ P '>=XT(p)xr(p')(^) d f e iL ^'> x dx(7) 

where \A is defined for any set A as \a (x) = 1 if x € A 
and xa(x) — otherwise. For the continuous case the 
mentioned inequality 0] is not helpful, since the Hilbcrt 
space associated with any set Q is infinite dimensional, 
so Tr e = oo . We proceed as follows: write instead 

h(t) < ey/t{l -t)- Ct(l - t) log e (8) 

valid for small enough e (with a different constant C). 
We then take trace of both sides. It remains to es- 
timate Tr y/PQP{\ - PQP). We have the following: 
Tr y/PQP(l - PQP) < Tr y/PQ a P where Q a is a pro- 
jection on a box containing Q, and we have used op- 
erator monotonicity of t — > i 1 / 2 (see e.g. ; 19|), and that 
PQ^P > PQP (as operators). Next we note that the op- 
erators PQ^P and QaPQa have the same positive eigen- 
values counted with multiplicities j^H. Thus we have 
Tr^PQP(l-PQP) < Tr VQnPQn < TrVQnPoQn 
where we have used the monotonicity again. It remains 
to evaluate Tr \/ QaPaQa this can be done using bounds 
on the singular values of the operator QrjfbQrj, which 
in this case are also the respective eigenvalues. It fol- 
lows from |23| that for any a > d/2 the eigenvalues 
of Q D P D Q a satisfy A » < CrT x l' 1 - a l d L d l' l+a . Taking 
a > 3d/2 we find Tr y/Q n PaQa < L d+S for any 6 > 0, 
and thus we can choose e.g. e < L~ A_1 , and © follows. 

Having established (AN) 2 as a way of obtaining 
bounds we proceed to evaluate it. Our next result is 

Theorem. Let f2,T be two compact sets in M. d , d > 1, 
with smooth boundaries dil, dT. Then 

+ o( J L d - 1 log J L) 

where n x ,n p are unit normals to <9f2, <9r, respectively. 
The full proof is too technical to be included here, and 
will appear elsewhere (see however |2l|). It starts by 
observing that 

z r \ 2d r 

Tr(PQP) 2 =( — ) / A n (z)A T (Lz)dz (10) 

V 27T/ jRd 

where An(#) = / Rd Xn( x )xn( x - z)dx is the volume of 
the set f2 intersected with shifted by z (i.e. f2n(S! + z)), 
and proceeds with an asymptotic analysis of this integral. 

Note the geometric nature of the coefficient in © : for 
a spherical Fermi sea T, and a convex region O this co- 
efficient is just the average cross-section of f2 over all 
directions. In the more general case the coefficient de- 
pends on the two surfaces and their mutual orientations. 
Thus © with © establish the scaling 0(L d ~ x logi) < 
S < 0(L d_1 (logL) 2 ) when T, fi have smooth boundaries. 



Fractal Boundaries.- A very interesting enhancement 
in the scaling of S occurs if the sets f2,r are allowed to 
have fractal-like boundaries. Physically, this means that 
making the boundaries less regular, makes the typical 
momentum states more incompatible with the shape of 
the region f2, and hence contribute to enhanced entropy 
when integrating the external modes. More precisely, it 
was shown in |l[ that if < Vol(ft \ (fi + h)) < 

C 2 \\h\\^ for small \\h\\ and some < Pn,p T < 1, and 
the same holds for T with p T , then (AN) 2 is bounded 
above and below by (7i, 2 £ d ~ min(/3f! ' /3r - ) if Pn + Pr and 
Ci^L^ 13 " logL if Pq = f3 T . In particular this and © 
imply that S > 4CiL d - & logi if /3 n = Pr and S > 
4C , ii d - min ( /3 "'' 9r ) if Pn ^ p r 

Connection to Widom's conjecture. - It turns out that 
the result © is a special case of a well-known conjecture 
by H. Widom (1982). The problem of time-frequency 
limiting mentioned in the introduction leads to a study 
of the spectrum of the operator PQP where Q is a time 
window scaled by L, and P represents a frequency win- 
dow. One way of studying the eigenvalues of PQP is to 
study the asymptotic behavior of Tr / (PQP), as L — > oo, 
for some general class of /. It is conjectured in 0] that 
for a function f(t), analytic on a disc of raduis > 1 with 
f(0) — 0, the following holds as L — > oo 

Trf(PQP)= (J-^ f(l) J J dxdp 

+ o(L d - 1 \ogL) 

where n x ,n p are unit normals to dVL, dT, respectively, 
and U(f) = /J dt. The formula (TTJ and a gen- 

eralized form of it were proved for d = 1 i n fl/^ an d [Til . 
For d > 2 only special cases were proved |2l|, |23 • Note 
finally that © is a verification of Widom's conjecture for 
the special case f(t) — t(l — t). 

In a broader context one may think of Widom's con- 
jecture (|llfl as a generalization of the strong (two-term) 
Szego limit theorem (SSLT) for the continuous setting. 
The SSLT plays a special role in entanglement entropy 
H3| The SSLT was initially used by Onsager in his 
celebrated computation of the spontaneous magnetiza- 
tion for the 2D Ising model (see e.g. [3). It is inter- 
esting to note that in Onsager's computation (and also 
in |2J]) the leading asymptotic term vanishes, and one 
needs to compute the sub-leading term. This is exactly 
the situation that we have in the continuous version of 
the Szego theorem l|ll[) : the leading term should vanish 
since h(l) = 0. 

Widom's conjecture suggests the explicit geometric ex- 
pression for the entropy QJ. Note that if £1 = T = [0, l] d 
then the double integral in 0) equals 4d (twice the num- 
ber of faces), so that © and © are consistent. Note also 
that the coefficient ^| in the expression for the number 
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variance © gives a lower estimate for S in within 
16% of the conjectured in QJ. 

Finite temperature.- From the semiclassical point of 
view one expects the entropy to be extensive, 5 ~ L d , 
for T > 0. This suggests to look for a transition tem- 
perature between the L d and Lr -1 log L regimes. Let 
[3 = l/T and introduce the Fermi-Dirac function k(p) — 
l/(l +e /3(bl 2 -M)) (we take ft = fc B = land m = 1/2). The 
expression for the entanglement entropy at finite temper- 
atures || is given by Tr h(QKQ) where K is the operator 
of multiplication by k(p) in momentum space. Semiclas- 
sically, integrating over the phase space one finds [sij l 



Tr h(QKQ) 



V27T 



Vol(O) / h(k(p)) dp + 0(L 



d-l\ 



Introducing polar variables and scaling out (3 gives 



x h 



1 



1 + u 



1 



\ogu\- 1 + d / 2 



(12) 



+ 0(L d - 1 ) 



which scales as L d p^ 1 fi~ 1+d ' 2 for (3 — ► oo. Comparing 
this with the T = results above we see that for the zero 
temperature effect L d_1 logL to be seen, the transition 
temperature should satisfy T^ 1+d / 2 ~ -^J— , L — > oo. 

Summary and Discussion. - In systems with finite cor- 
relation length £, one expects quantities such as the en- 
tropy S and the number variance (AiV) 2 to scale like the 
area of the boundary of the region. The system stud- 
ied here does not behave this way. Here, the correlation 
function < a\.a x i > decays slowly and the fermion mo- 
mentum modes are spread over the entire system and are 
highly sensitive to localization in space and consequently 
the area law is violated. 

Let us summarize the concrete results of this Letter: 
We prove that the scaling is of the form L d ~ 1 \ogL for 
cube like domains. We find a connection between the 
scaling behavior of S and a well-known conjecture due to 
Widom i|llfl . which suggests the explicit geometric for- 
mula Q for S in any d. Finally, while Widom's conjec- 
ture is far from being proven, we find that it holds for 
(A7V) 2 , and use this to obtain lower and upper bounds 
on S. We also find an enhanced scaling of S for fractal 
like boundaries and at finite temperatures. 
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We believe the last argument can be made rigorous, but 
that is outside the scope of the present paper. The func- 



tion 
of 



-t) 



, , — — is not Riemann integrable and hence the result 
(l7| is not applicable as such. 
Note that JJJ coincides with the upper bound in JjjJ. 
In order to justify the interchange of the asymptotic limit 
in L, and the series expansion. 

This holds since for any positive trace class operators 
A, B the fact that Tr A" = TrB™ for all n = 1,2, ••■ 
implies that A and B have the same positive eigenvalues. 
A fractal set similar to the one appearing in [2l| was 
independently constructed for d — 1 in fl^l . 
Indeed, for translation invariant systems, the correlation 
matrices are Toeplitz matrices (i.e. < ajoj >= g(i — j)), 
and the asymptotics of Toeplitz determinants are given 
by (various versions of) the SSLT. 

In |2a| a three-term asymptotics of Tr f(QKQ) is proven 
for analytic J. However 1121 strictly speaking does not 
follow from 1251 since h is not analytic. 



